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We investigate the phase diagram and the spin-orbital entanglement of a one-dimensional 
SU(2)®XXZ model with SU(2) spin exchange and anisotropic XXZ orbital exchange interac¬ 
tions and negative exchange coupling. As a unique feature, the spin-orbital entanglement entropy 
in the entangled ground states increases here linearly with system size. In the case of Ising or¬ 
bital interactions we identify an emergent phase with long-range spin-singlet dimer correlations 
triggered by a quadrupling of correlations in the orbital sector. The peculiar translational invariant 
spin-singlet dimer phase has finite von Neumann entanglement entropy and survives when orbital 
quantum fluctuations are included. It even persists in the isotropic SU(2)0SU(2) limit. Surpris¬ 
ingly, for finite transverse orbital coupling the long-range spin singlet correlations also coexist in 
the antiferromagnetic spin and alternating orbital phase making this phase also unconventional. 
Moreover we also find a complementary orbital singlet phase that exists in the isotropic case but 
does not extend to the Ising limit. The nature of entanglement appears essentially different from 
that found in the frequently discussed model with positive coupling. Furthermore we investigate 
the collective spin and orbital wave excitations of the disentangled ferromagnetic-spin/ferro-orbital 
ground state and explore the continuum of spin-orbital excitations. Interestingly one finds among 
the latter excitations two modes of exciton bound states. Their spin-orbital correlations differ from 
the remaining continuum states and exhibit logarithmic scaling of the von Neumann entropy with 
increasing system size. We demonstrate that spin-orbital excitons can be experimentally explored 
using resonant inelastic x-ray scattering, where the strongly entangled exciton states can be easily 
distinguished from the spin-orbital continuum. 

PACS numbers: 75.25.Dk, 03.67.Mn, 05.30.Rt, 75.10.Jm 


I. INTRODUCTION 

Spin-orbital coupling phenomena are ubiquitous in 
solids and have been known to exist since the early days 
of quantum mechanics and band theory, but only re¬ 
cently it was realized that the quantum nature of or¬ 
bital degrees of freedom plays a crucial role in the fields 
of strongly correlated electrons BH3 and cold atoms @- 
[l2:]. The growing evidence of spin-orbital entanglement 
(SOE) accumulated due to novel experimental techniques 
which probe a variety of underlying electronic states. The 
strong Coulomb interactions and the relativistic spin- 
orbit interaction entangle locally the spin and orbital de¬ 
grees of freedom [13] which display an amazing variety 
of fundamentally new and fascinating phenomena, rang¬ 
ing from topologically nontrivial states 0 , relativistic 
Mott-insulating behavior in 5 d [HI, [HI and 4 d mm 
transition-metal oxides and entanglement on superex¬ 
change bonds in spin-orbital models @,0. Other more 
recent developments include entangled spin-orbital exci¬ 
tations Um doped spin-orbital systems [22j , skyrmion 
lattices in the chiral metal MnSi (23|, multiferroics, spin- 
Hall effects [24], Majorana and Weyl fermions (25|, topo¬ 
logical surface states [26j , Kondo systems [27j , exotic spin 
textures in disordered systems, to name just a few. 

To date, experimental observation of a dynamic spin- 
orbital state has been a challenge. Apart from the in¬ 


trinsic anisotropy and the relative complexity of the or¬ 
bital couplings, it has been shown that the interplay be¬ 
tween the two frustrated degrees of freedom may lead 
to exotic states of matter. An x-ray scattering study 
of a dynamic spin-orbital state in the frustrated magnet 
Ba 3 CuSb 2 0g supports spin liquid state [H, [29], while 
FeSc 2 S 4 |30l - l32l and the d 1 effective models on the trian¬ 
gular lattice [33j] and on the honeycomb lattice ii 
are found to be candidates for spin-orbital liquids in the 
theory. Recently remarkable progress was achieved due 
to rapidly developed resonant inelastic x-ray scattering 
(RIXS) techniques [HI which helped to explore the el¬ 
ementary excitations in S^CuOs m Hi and Sr 2 lr 04 
[39| . with antiferromagnetic (AF) and ferro-orbital (FO) 
order in ground states. Orbital order in the spin-gapped 
dimerised system SrsC^Og below the Jahn-Teller tran¬ 
sition was also identified [HI- However, it remains chal¬ 
lenging experimentally and theoretically, mainly owing 
to the lack of an ultimate understanding of spin-orbital 
correlations. 

In the Mott insulators with an idealized perovskite 
structure, the low-energy physics is described by spin- 
orbital models, similar to the Kugel-Khomskii model [2], 
where the spin and orbital are considered on equal foot¬ 
ing as dynamic quantum variables Q. Spin interaction 
possesses SU(2) symmetry, which will be broken however 
by the relativistic spin-orbit coupling. It couples spins 
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to the orbitals, that are in general non-SU(2)-symmetric 
in a solid. However, this coupling can frequently be ne¬ 
glected in realistic 3d systems and one is left in general 
with entangled spin-orbital superexchange problem [gj, 
that is the eigenstates cannot be written as products of 
spin and orbital wave functions. One immediate conse¬ 
quence of entanglement is that spin and orbital terms 
cannot be factorized in the mean-field approach. Or¬ 
bitals are spatially anisotropic and thus their interactions 
have lower symmetry than the spin ones which reflects 
the directional dependence of the orbital wave functions. 
For the fixed occupation of orbitals, the magnitude and 
sign of the spin-orbital superexchange interactions follow 
the classical Goodenough-Kanamori rules [4l[, but quan¬ 
tum fluctuations change them and make it necessary to 
consider spin-orbital interplay in entangled states on ex¬ 
change bonds 0- Therefore, it is important to measure 
whether eigenstates are entangled or not. 

A natural measure of SOE is the von Neumann entropy 
(vNE) which we write first for the nondegenerate ground 
state |T 0 ), 

SvN = -Tr A {^ 0) log 2 ^ 0) }. (1.1) 

Here we consider a system D composed of two non¬ 
overlapping subsystems |42|, i.e., H = dUfi, dflfi = 0, 
and is the reduced density matrix. It is obtained 
by integrating the density matrix over subsystem H, i.e., 
pfSj) _ Ttb|To)(To|. However, one has to realize that 
information contained in entanglement entropy depends 
crucially on how one partitions the Hilbert space of the 
system. To investigate SOE we use here as two subsys¬ 
tems A and B the spin and orbital degrees of freedom in 
the entire chain. Standard spin-orbital phases may have 
entanglement in only one sector and here we concentrate 
on joint SOE [43j. This choice is distinct from the one 
conventionally made when the system is separated into 
two spatially complementary parts Q, for instance in 
frustrated spin chains [45| or in the periodic ID Ander¬ 
son model J46j. 

Though much attention was devoted to the ground 
state in the past E3, it has been noticed only recently 
that the entanglement entropy of low-energy excitations 
may provide even more valuable insights [43|, HH, Hgj 
which are of crucial importance to understand the origin 
of quantum phase transitions in spin-orbital systems [50j . 
The well known area law of the bipartite entanglement 
entropy restricts the Hilbert s pac e accessible to a ground 
state of gapped systems (fill, [52j, while the area law is 
violated by a leading logarithmic correction in critical 
systems, whose prefactor is determined by the number of 
chiral modes and precisely given by Widom conjecture 
[53| . In this respect, the application of the entanglement 
entropy in describing quantum criticality in many-body 
Hamiltonian merits a lot of studies [ 42 , [54| . 

On the other hand, the excited states have the mix¬ 
ture of logarithmic and extensive entanglement entropy, 
and the logarithmic states are expected to be negligible 


in number compared to all the others. The entanglement 
in excited state is proven always larger than that of the 
ground state of a spin chain j55|. For a spin-orbital cou¬ 
pled system, the division of spin and orbital operators 
retains the real-space symmetries, which is beneficial to 
the calculation of mutual entanglement. In two-particle 
states, the SOE is determined by the inter-component 
coherence length [43j| , as though the state has sufficient 
decay of correlations [45| . 

The aim of this paper is to use the entanglement en¬ 
tropy to investigate the full phase diagram of the one¬ 
dimensional (ID) anisotropic spin-orbital SU(2 )<g>XXZ 
model. The main motivation for considering the Ising 
asymmetry in the orbital sector comes from the obser¬ 
vation that spin-orbital entanglement is large when both 
subsystems, i.e., spin and orbital sectors, reveal strong 
quantum fluctuations. Thus the Ising anisotropy which 
is present in many physical systems introduces addi¬ 
tional control of orbital fluctuations and thereby pro¬ 
vides an important control parameter for SOE. Here 
we focus on the model with negative exchange inter¬ 
action. This choice of the exchange coupling restricts 
somewhat joint spin-orbital fluctuations being particu¬ 
larly large near the SU(4) symmetric point in the ID 
spin-orbital model with positive coupling constant [ 56 | . 
but opens novel possibilities for entangled states, as we 
show below [50]. An interesting phase with entangled 
ground state, consisting of alternating spin singlets along 
the spin-orbital ring, is found for Ising orbital interac¬ 
tions when the dimerization in the spin channel induces 
the change from FO to alternating orbital (AO) correla¬ 
tions. Here we report the complete phase diagram of the 
anisotropic SU(2 )<g>XXZ spin-orbital model, with two 
phases of similar nature which gain energy from singlet 
correlations leading to dimerization, either in spin or in 
orbital sector. These phases were overlooked before in 
the fully symmetric case, i.e., in the phase diagram of 
the isotropic SU(2)®SU(2) model (43j . 

We also analyze the nature of spin-orbital excited 
states, particularly in the case of the disentangled ferro¬ 
magnetic (FM) and FO ground state, labeled as FM/FO 
order. We also analyze entanglement entropy in the ex¬ 
cited states for the FM/FO phase and show that spin- 
orbital excitations form a continuum, supplemented by 
collective bound states. The latter states are character¬ 
ized by a logarithmic scaling behavior, and as we show 
could be detected by properly designed RIXS experi¬ 
ments (EM- 

The paper is organized as follows. The model is intro¬ 
duced in Sec. HD In Sec. inn we present an analytic solu¬ 
tion for the ground state in the Ising limit of the orbital 
interactions. A more general situation with anisotropic 
XXZ orbital interaction is analyzed in Sec. ESI and 
the phase diagram for the isotropic SU(2)®SU(2) model 
is reported in Sec. IIVB1 This model and the obtained 
SOE are different from the AF case, as shown in Sec. 
II V Cl Next we determine the elementary excitations in 
the FM/FO phase in Sec. IVIand show that they are en- 
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tangled although the ground state is disentangled. The 
vNE spectral function is presented in Sec. IVI Ai includ¬ 
ing the scaling behavior of the bound states which is 
contrasted with that in the AF/AO ground state. In 
Sec. IVI Bl we explore the possibilities of investigating 
entanglement in the present ID spin-orbital model by 
RIXS. The paper is concluded by a discussion and brief 
summary in Sec. IVII1 Some additional technical insights 
which are accessible by an exact solution of the two-site 
model are presented in the Appendix. 

II. THE ID SPIN-ORBITAL SU(2 )®XXZ 

We consider the ID spin-orbital Hamiltonian which 
couples S = 1/2 spins and T = 1/2 orbital (pseudospin) 
operators, 

A,y), (2.1) 
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with SU(2) spin Heisenberg interaction 'Hj(x), orbital 
anisotropic XXZ interaction Hj(A, y), 

H S j ( x ) = Sj • Sj+ 1 + x, (2.2) 

nj (A ,y) = A (T/T/ +1 + TJTJ +1 ) + T/T/ +1 + y. (2.3) 

We take below J = 1 as the energy unit. The model 
Eq. (j2.1| ) has the following parameters: (i) x and y 
which determine the amplitudes of orbital and spin ferro- 
exchange interactions, — Jx and — Jy , respectively, and 
(ii) A which interpolates between the Heisenberg (A = 1) 
and Ising (A = 0) limit for orbital interactions. When 
A = 1, the spin and orbital interactions are on equal 
footing and the symmetry of the Hamiltonian (|2.1|) is en¬ 
hanced to SU(2)0SU(2) — this model describes a generic 
competition between FM and AF spin, and between FO 
and AO bond correlations [43[. 

We emphasize that the coupling constant — J is neg¬ 
ative , so at large x > 0 and y > 0 it gives a disentan¬ 
gled FM/FO ground state, see below — therefore the 
model may be called in short FM. This choice of the 
exchange coupling restricts somewhat joint spin-orbital 
fluctuations being large near the SU(4) symmetric point, 
(x,y) = (0.25,0.25), in the ID spin-orbital model with 
positive , i.e., AF coupling constant (56|, but opens other 
interesting possibilities for entangled states, as we have 
shown recently [5(J. Both total spin magnetization S z 
and orbital polarization T z are conserved, and time rever¬ 
sal symmetry leads to the total momentum either k = 0 
or k = 7r. 

Before analyzing the spin-orbital model of Eq. (1) in 
more detail, let us summarize briefly the properties of 
the well known AF model, with positive coupling con¬ 
stant J. The SU(4) symmetric Hamiltonian found at 
(x,y) = (0.25,0.25) is an integrable model which can 
be solved in terms of the Bethe Ansatz m 6oj. Away 
from the SU(4) symmetric points this choice of the cou¬ 
pling constant favors the phases with spin-orbital order 


depending on the actual values of x and y, and the phase 
diagram obtained by numerical methods includes in gen¬ 
eral phases with all types of coupled spin-orbital order, 
i.e., FM/FO, AF/FO, AF/AO, and FM/AO, as well as 
the gapless spin-orbital liquid phase near the SU(4) point 
[6li l62|. In addition, Schwinger boson analysis gives 
phases with spin-orbital valence-bond correlations and 
also spin valence bond and orbital valence-bond phases 
(63|. The latter two show a tendency towards dimerised 
spin or orbital correlations which occur here in the prox¬ 
imity of the SU(4) point. For some special choice of pa¬ 
rameters the model can be solved exactly: (i) when A = 1 
and x = y = 3/4, the exact ground state is doubly de¬ 
generate with the spins and the orbitals forming singlets 
on alternate bonds, while (ii) when A = 0, x = 3/4 and 
y = 1/2, the non-Haldane spin-liquid ground state can 
be analytically obtained [63, |65j , and (iii) several inte¬ 
grable cases were presented for interactions with special 
symmetries Km 1 6 71] , or (iv) with XY orbital interactions 

(a = oo) Ear 

The form of Eq. (|2.ip is not the most general one 
but is representative for real spin-orbital systems with 
anisotropic orbital interactions. In real systems the or¬ 
bital part contributes by additional superexchange terms 
which are not coupled to SU(2) spin interaction 14]. For 
instance, in the case of t^g orbital degrees of freedom as 
in the perovskite titanates or vanadates, the interactions 
along the c cubic axis involve the doublet of two orbitals 
active along it, i.e., the yz and zx orbitals [Hj; a similar 
situation is encountered in a tetragonal crystal field of a 
quasi-ID Mott insulator [6jJ, or for p x and p y orbitals of 
a ID fermionic optical lattice HClq|. 

A priori , due to the quartic spin-orbital joint term, 
(X (yy.iiXfYjvy., -Ypyi-'/JY;.,] in the Hamil¬ 
tonian Eq. m the spin-orbital interactions are entan¬ 
gled, and the spin and orbital operators cannot be sepa¬ 
rated from each other in the correlation function, except 
for some ground or excited states in which the SOE van¬ 
ishes. The spin-orbital bond correlations (J2.4D 

cr * ((Sj ■ S j+ t )[A (TjTj +1 + T/Tjq) + r; T? +1 ]) , 

(2.4) 

are uniform in the considered system and C\ ot does not 
depend on the site index j. We investigate below these 
composite quartic correlations and show that they could 
also be surprisingly large. As an additional criterion of 
setting up the phase diagram, we use below the fidelity 
susceptibility which elucidates the change rate of ground 
states in the parameter space |J0]. It serves as an or¬ 
der parameter to characterize the phase diagram of the 
anisotropic (A < 1) spin-orbital model (12.ip . The fidelity 
susceptibility is defined as follows, 

/w O l- ln.F(A,(5A) 
w (A) = -2hm (sx)1 , (2.5) 

where the fidelity 

7-(A,«A) = |($o(A)|®o(A + JA))|, 


( 2 . 6 ) 
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is taken along a certain path in the parameter space in 
the vicinity of the point A = A(A ,x,y). 


III. ISING ORBITAL INTERACTIONS (A = 0) 

In the Ising limit of orbital interactions (A = 0) the 
Hamiltonian (|2.1j) simplifies and has SU(2)(g)Z2 symme¬ 
try — it is a prototype model for the directional orbital 
interactions with quenched quantum fluctuations in t^g 
systems. This may happen in real compounds in two 
ways: (i) either only one of the two active orbitals is oc¬ 
cupied by one electron and contributes in the hopping 
processes along the 180° bonds [zH or 90° bonds [72|, 
or (ii) the orbital degrees of freedom are quenched in 
the presence of strong crystal field. In both these cases 
the orbital exchange (orbital-flip) processes are blocked 
and orbital interaction are of a classical Ising-like form. 
Such Ising interactions are frustrated when they emerge 
in higher dimension, as in the well-studied orbital com¬ 
pass model [731475) and in Kitaev model [76], see also a 
recent review on the compass model It is now in¬ 

triguing to ask what happens to the SOE in this case. 
It may be still triggered by spin fluctuations while the 
model with Ising spin interactions (fpi) is classical. 

The phase diagram of the model Eq. (G3D at A = 0, 
i.e., in the absence of orbital fluctuations, which follows 
from fidelity susceptibility (|2.5|) is displayed in Fig. [lj 
As expected, one finds four trivial combinations of spin- 
orbital order: FM/FO (phase I), AF/FO (phase II), 
AF/AO (phase III), and FM/AO (phase IV). All these 
phases have the entanglement entropy (HD 

^vN - 0 

and spins and orbitals disentangle. Transitions between 
pairs of them are given by straight lines and may be also 
obtained rigorously by the mean-field approach. The 
ground state of a L-site chain stays in the subspace 
S z = 0 , T z = 0 , momentum k = 0 (always degenerate 
with S z = 0, T z = 0, k = 7 r for all parameters) in phases 
III (AF/AO), IV (FM/AO) and V, while it is found in 
the subspaces S z = 0, T z = ±L/ 2 , k = 0 in phases I 
(FM/FO) and II (AF/FO) (of course, in phases I and IV 
also other values of S z 7 ^ 0 , with — L/2 > S z > L/ 2 , are 
allowed and the ground states have the respective degen¬ 
eracy) . The ground states with energy Eq = 0 are highly 
degenerate when x < —1/4 along the critical line y = 1/4 
between phases III and IV, suggesting that antiparallel 
orbitals erase the spin dynamics. Along the critical line 
y = —1/4 between phases I and II, the ground states are 
also highly degenerate when x > 3/4, and parallel or¬ 
bitals on the bonds (in FO order) quench again the spin 
fluctuations. 

Although the orbital interactions are Ising-like, entan¬ 
gled spin-orbital ground state occurs in phase V. In order 
to understand better emergent phase V, we introduce the 
longitudinal equal-time spin/orbital structure factor, de¬ 
fined for a ring of length L (with a lattice constant a — 1; 



FIG. 1. (Color online) Spin-orbital entanglement entropy <S° N 
Eq. CH) and the phase diagram in the ( x , y) plane of the 
SU(2)0Z2 spin-orbital model (12.11) with A = 0 as obtained for 
the system size of L = 8 sites. The critical lines are discerned 
by both fidelity susceptibility and analytical method. Phases 
I-IV are disentangled (<S° N = 0) with order defined as follows: 
FM/FO (phase I), AF/FO (phase II), AF/AO (phase III), and 
FM/AO (phase IV). The spin and orbital textures in phase 
V with finite entropy <S° N > 0 are explained in the text. 


we use periodic boundary conditions) by 

1 L 

S zz {k) = - e~ ik U-i\S z Sp), (3.1) 

1 L 

T zz {k) = — \TjTp). (3.2) 

jJ'= 1 

The calculation of the equal-time structure factor S zz (k ) 
for a model of uncorrelated nearest neighbor dimers 
was compared with the one for the kagome lattice 
ZnCu 3 (OD) 6 Cl 2 [78j]. One finds analytically that in the 
case A = 0, a cosine-like spin structure factor, i.e., 
S zz (k ) oc (1 — cos &), is revealed in phase V for y = —1/4, 
implying that only nearest neighbor spins are correlated. 
This finding is essential as the short-range spin correla¬ 
tion indicates here a translation invariant dimerised spin- 
singlet state which has the same spin structure as the 
Majumdar-Ghosh (MG) spin state [79). However, this 
state is not triggered here by frustrated interactions J\ 
and J 2 , but is evidently induced by the correlations in 
the orbital sector. 

In the Ising limit we obtain the analytic ground state 
for phase V as described below. The essential feature is 
that the energy is gained by spin singlets occupying the 
bonds with AO states, while the bonds connecting two 
spin singlets have FO order, see Fig. [2ja). To construct 
the ground state, we introduce the corresponding four 
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(a) 

t *Ht t» ♦ * 


(b) 

f +► t<+ t +►<♦ * 


(c) 


FIG. 2. (Color online) (a) One of four translational equiva¬ 
lently spin and orbital configurations in the Ising limit of the 
spin-orbital model (12.11) at A = 0 and y — —0.25. The spins 
form isolated dimers (shaded ovals), (b) A single orbital ex¬ 
citation and induced spin configuration, (c) A single spin flip 
makes a singlet-triplet spin excitation, but does not induce 
any change in orbital correlations. 


+ _ .+ + _ .+ 
M> M t * ♦ ♦ 



+_ .++. _+ 


configurations in the orbital sector: 

|0l) — | + d-M-), 

102) = 1- + + -- + +—), 

| 0 3 ) = | 

104) = | H-H-1-)• (3-3) 

The solutions are classified by the momenta correspond¬ 
ing to the translational symmetry of the system. The 
orbital wave functions in the ground state for the mo¬ 
menta k = 0, 7t/2, 37t/2, 7t correspond to: 

10*) = \ (|0i) + e ife |0 2 > + e 2ife |0 3 > + e 3ifc |0 4 » • (3.4) 

In spin subspace there are two distinct (but nonorthogo- 
nal) states: 


The state |^/c=o) is entangled both in individual spin and 
orbital subspaces, and also is characterized by SOE along 
the chain. Such a many-body state, and similar states ob¬ 
tained for other momenta, k = ±ir/2 and k = 7r, give an 
exact value of the vNE, = 1. The resulting fluctu¬ 
ations between these states suppress conventional order 
and the system features finite entropy even at zero tem¬ 
perature, in contrast to the naive expectation from the 
third law of thermodynamics. The emergent excitations 
are also entangled and fundamentally different from the 
individual spin or orbital ones, see Figs. [2jb) and[2fc). 

As the orbital correlations are classical in the Ising 
limit, we can determine all the phase boundaries analyt¬ 
ically by considering the spin interactions for various or¬ 
bital configurations. The lower boundary between phase 
III (AF/AO) and V at 1 / = -1/4 can be determined by 
comparing the uniform state with energies of AO corre¬ 
lation on a bond, i.e., (T?T? +1 ) = —1/4, with the al¬ 
ternating state of pairs of the same orbitals shown in 
Fig. [2fa), i.e., (T^T^ +1 ) = (—l) J /4, which coexists with 
spin dimer order (spin interactions vanish for a pair of 
identical orbitals). One finds the following effective spin 
Hamiltonian in this case: 

H d1 m = \J2 (3.7) 

j'Godd 

and the corresponding ground state energy per site in the 
thermodynamic limit is 

^dim = ^ (—0.75 + x ). (3.8) 

The dimerised phase competes with the AO order coex¬ 
isting with the ID resonating valence-bond spin state, 
with energy 


£ao = 1 (-0.4431 + ®). (3.9) 

Hence, one finds that IM < E^ Q for x > 0.136. The 
quadrupling due to spin-orbital interplay in phase V is 
well seen by the calculation of the four-spin correlation 
function which we define following Refs. |80l. l81j. 


iV'f) = [1,2][3,4] - - - [AT - 1, AT], 

ItT) = [2,3][4,5] • • • [V, 1], (3.5) 

where the singlets are located on odd (even) bonds. Here 
a singlet is defined by [/, / + 1] = (| ti) — | 

One representative component of the ground state with 
the orbital part \(j)^) accompanied by the spin state \^i) 
is shown in Fig. [2fa). The ground state in the k = 0 
subspace is given by the superposition 


D(r) = [((£ ■ S i+1 )(S i+r • S i+r+ 1 )) 

i 

If y = —1/4, spin dimer correlations alternate and 

2 


(3.10) 


D(r) = (-1 y 


(3.11) 


|$k= 0 ) = 


2 (I0l)®|02 > ) + 102)0 |0?) + |03><S ) |'02 3 ) + 104)® |0?)) 

= T (M±JM 0 \^) + M±JM 0 tyD)) .(3.6) 


y/2 V y/2 


V2 


which follows from Eq. (|3.10[) for the alternating spin 
singlets, (Si • Si+ r ) = — 3[1 — (—l) r ]/8. Indeed, one finds 
this value (|3.11| ) for x G [0.2,0.7] and the result is robust 
and the same for systems sizes L = 12 and L = 16, see 
Fig. [3j On the contrary, for x < 0.2 the values of D(r) 
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FIG. 3. Dimer correlation function D(r) (13.101) obtained for 
the anisotropic SU( 2 )®Z 2 spin-orbital model for different val¬ 
ues of x £ [0, 0.7] in phases V and III and for the ring of length: 
(a) L m 12, and (b) L — 16 sites. Parameters: y = —0.25 and 
A = 0. 


decrease with increasing distance r, and would vanish in 
the thermodynamic limit of L —>■ oo. 

When y < —1/4, there are three competing phases 
with predetermined orbital configurations (AO, DIM, or 
FO) and the corresponding spin interactions given by ef¬ 
fective spin Hamiltonians: 

ffAO =(\~y) £($• • Sj+ 1 + X), (3.12) 

' ' 3 



- (j + v\ £ (Sj • Sj+ 1 + x), (3.13) 

' ' j'Geven 

H FO = -(\+y) £(5,- • S m + x). (3.14) 

' ' 3 


In this case also even inter-singlet bonds contribute to the 
energy in the |DIM) state, but the spin correlations van¬ 
ish, i.e., (Sj • Sj+ 1 ) = 0. It is obvious that Hao and F7 fo 
stand for the same (translational invariant) spin Hamilto¬ 
nian, and Ufo will have lower ground state energy when 
x > —(Sj • Sj+ i)af — 0.4431. The dimerised AF Heisen¬ 
berg chain (|3.13j) related to spin-Peierls state cannot be 
solved trivially, with the exception of the free-dimer limit 
(y = —0.25) and the uniform Heisenberg limit (y = —oo) 
[82|, One finds the ground state energy per site £oc(d) of 
a pure dimerised spin chain (83|, 


£oo(£) 


3 1 

4 1 + a 


a 2 a 3 

1 T — T — T ■ 

8 32 


(3.15) 


with a = (1 — (5)/(l + <5) and 5 = l/|4y| > 0.4. For 
<5 < 0.4, 

e 00 (^)=ln2-(ln2-l)|(5| 4 / 3 . (3.16) 


In such a case, E^ lu = y[s 00 (S) — x\. The overwhelming 
dimerised phase will persist in a range of negative values 
of y , and the boundaries close at y = — oo, as is indicated 
by structure factors and fidelity susceptibility. 

The phase transitions in the phase diagram of Fig. [1] 
imply the discontinuous changes of order parameters in 
first-order quantum phase transitions. The orbital order 
changes from phase II (AF/FO) to phase III (AF/AO), 
as shown in Ref. [Ho|, but the Neel order persists in both 
of them and manifests itself in the two-spin correlation, 
(SfSi +r ). For translational invariant and orthonormal 
linear combinations of the symmetry-broken Neel (AF) 
states, 

i$y> = i nn • • • n>, 

l$2 F > = l (3.17) 

there are spin (SfSf +r ) = (—l) r /4 and dimer D(r) = 0 
correlations (for r / 0), while for dimer states |4>P IM ) 
and |<£> 2 ™), (SfSf +r )=0 (for r/±1) and D(r) ^ 0 (for 
r / 0), see Fig. [SJ respectively. These results reflect 
the long-range nature of the two types of order. The AF 
classical spin correlations are replaced by a power 

law for the AF spin 5 = 1/2 chain in the thermodynamic 
limit, 


(Si ■ s i+r ) ~ (-ir^H. (3.18) 

which is equivalently revealed by the structure factors 
S zz (k ) and T zz (k ) defined by Eqs. (|3.1| ) and (|3.2[) . 


IV. ENTANGLEMENT IN THE GROUND 
STATES 

A. The anisotropic orbital interactions (0 < A < 1) 

When A = 0, there is no dynamics in the orbital sec¬ 
tor, and the orbital structure factor is dominated by a 
single mode which follows from the Z 2 symmetry [50]. 
This changes when 0 < A < 1 and the quantum fluctu¬ 
ations in the orbital sector contribute. In order to un¬ 
derstand the modifications of the phase diagram in the 
entire interval 0 < A < 1 , we select A = 0.5 and study 
the longitudinal equal-time spin/orbital structure factor, 
defined for a ring of length L in Eqs. (|3.ip and (|3.2p . 

The most important change at finite A occurs for the 
phase transition between phases V and III (AF/AO) 
which becomes continuous for fixed y and decreasing x, 
with a gradual change of spin correlations from the al¬ 
ternating singlets to an AF order along the chain [Hoj . 
Here we discuss in more detail the intermediate case of 
A = 0.5. First we address the phases with uniform spin- 
orbital order. The spin structure factors has distinct 
peaks at k = 0 for FM order and at k = tt for AF order. 
Similarly, one finds a maximum of the orbital structure 
factor T zz (k ) at k = 0 for FO order and at k = tt for AO 
order. These structure factors complement one another 
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and one finds that the spin correlations are somewhat 
weaker due to stronger spin fluctuations, while the or¬ 
bital fluctuations are moderate at this value of A = 0.5. 

The dimerised phase V is characterized by a remark¬ 
ably different behavior, see Figs. |U^c) andlU^d). Phase 
V, found at the E point in Fig. [5j has spin dimers ac¬ 
companied by an orbital pattern with the periodicity of 
four sites, see Fig. Ha). Spin correlations give a sharp 
maximum of S zz {k ) at k = tt as for AF states, while two 
symmetric peaks of T zz (k) at k = tt/2 and k = 37r/2 in¬ 
dicate quadrupling of the unit cell in the orbital channel. 
When the model evolves towards the SU(2)®SU(2) limit 
with increasing A, one expects also a similar phase VI 
with interchanged role of spin and orbital correlations. 
Indeed, this complementary phase emerges already at 
small A > 0 and is identified by the respective structure 
factors shown also in Figs. |U^c) and 0(d). 

The above analysis of the structure factors demon¬ 
strates that the phase diagram found for A = 0.5 con¬ 
tains six distinct phases, see Fig. [5] The quantum phase 
transitions at the boarder lines I-V and II-V are of first 
order. The phase transition between phases III (AF/AO) 
and V is a first order transition only for A = 0, and here 
this transition is continuous [Id]. As described above, 
phase VI emerges at finite A but is still quite narrow in 
the phase diagram of Fig. [5j Also the phase transition 
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2tc 
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FIG. 4. (Color online) The spin S zz (k ) (a,c), and orbital 
T zz (k ) (b,d) structure factors obtained for the selected points 
shown in the phase diagram of the anisotropic spin-orbital 
SU(2 )®XXZ model, see Fig. [5j (a,b) A, B , C and D in 
phases I-IV, and (c,d) E and F in phases V and VI. Param¬ 
eters: A = 0.5 and L — 8 sites. 


S vN 

0.0 



FIG. 5. (Color online) Phase diagram in the (x,y) plane and 
spin-orbital entanglement <S°n (jl-lj) (right scale) in different 
phases (indicated by color intensity) of the anisotropic spin- 
orbital SU(2 )®XXZ model (12.11) with A = 0.5, obtained for 
a ring of L — 8 sites. Phases I-IV correspond to FM/FO, 
AF/FO, AF/AO, FF/AO order in spin-orbital sectors. Or¬ 
bitals (spins) follow a quadrupled pattern accompanied by 
spin (orbital) dimer correlations in phase V (VI). Six labeled 
points are: A = (—1,1), B — (—1,-0.5), C = (0.5,1), 
D = (0.5,-0.5), E = (1, —0.5), and F = (-1,0.38) — they 
are used to investigate spin and orbital structure factors in 
different phases, see Fig. [4] The phase boundaries, deter¬ 
mined by the dominant modes of structure factors are shown 
by solid (dashed) lines for the first (second) order quantum 
phase transitions. 


from phase III to phase VI is continuous. We have ver¬ 
ified that due to the short-range nature of spin-orbital 
correlations, the size L = 12 is sufficient as the phase 
boundaries are here almost the same as for the ring of 
L = 8 sites. 


TABLE I. The spin-orbital configurations, momenta k , the 
total spin S z and orbital T z quantum numbers of the ground 
states I-VI found for the anisotropic SU(2 )(&XXZ model at 
A < 1. All states in these subspaces are nondegenerate (d — 
1), but in case of S z = L/2 there are L + 1 degenerate states 
for total S = L/2, and for T z — L/2 and A < 1 an equivalent 
state for — orbitals has T z — — L/2. At A = 0 phase VI is 
absent and the ground state degeneracy of phase V changes 
to d = 4 corresponding to momenta k — 0, ±7r/2, 7t. 


phase 

spin state 

orbital state 

k 

s z 

rj-iZ 

I 

tttttttt 

++++++++ 

0 

L/2 

L/2 

II 

nnnn 

++++++++ 

0 

0 

L/2 

III 

tmtm 

+ 
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+ 
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+ 
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+ 

1 

0 

0 

0 

IV 

tttttttt 

+ 

1 

+ 

1 

+ 

1 

+ 

1 

0 

L/2 

0 

V 

(S = 0) singlets 

+--++-+- 

0 

0 

0 

VI 

mtmt 

(T = 0) singlets 

0 

0 

0 
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The spin-orbital phases I-VI found at A > 0 are sum¬ 
marized in Table I. Phases I-IV have polarized or alter¬ 
nating spin and orbital components, combined in all pos¬ 
sible ways into phases: FM/FO, AF/FO, AF/AO, and 
FM/AO. Phases with either S z = L/2 or T z = L/2 
(FM or FO) have of course also degeneracy with respect 
to other possible values of S z or T z (the latter only at 
A = 1 when T is also a good quantum number). In addi¬ 
tion, there are two phases with dimer orbital (phase V) 
or dimer spin (phase VI) correlations. It is remarkable 
that these two phases survive in the isotropic model at 
A = 1, see Sec. HVBl 

For A > 0 also phase III is characterized by finite 
SOE, and it expands to higher values of y along vertical 
lines for fixed x < —0.25. Indeed, the onset of entangled 
region in the phase diagram of Fig. [5] moves to higher 
values of y with increasing A, see Fig. [6] At A = 0.25 the 
entanglement entropy S® N develops a narrow peak with 
a maximum at y ~ 0.15. This maximum broadens up 
and moves somewhat to the right (to higher y) when the 
orbital fluctuations increase with increasing A towards 
the isotropic SU(2)®SU(2) model. A sharp increase of 
the entropy to <S^ N = 1, visible in all the curves shown 
in Fig. [6j signals the SOE in phase VI which increases 
with increasing A. This large SOE can develop because 
phase VI is similar to phase V — it does not break trans¬ 
lation invariance of the model and different spin-orbital 
configurations contribute simultaneously. 

To detect SOE we employ here not only the vNE, 

Eq. (jl.ip . but also a direct measure by the spin-orbital 



y 


FIG. 6. (Color online) Spin-orbital entanglement entropy <S° N 
( 12 . 11 ) in the ground state of the spin-orbital model ( 12 . 11 ) as a 
function of y for selected values of A. The onset of phase 
VI is detected at A > 0 by a step-like increase of entropy to 
<S°n — 1- Parameters: x = —0.5 and L — 8. 



FIG. 7. (Color online) The onset of phase V and its gradual 
change into phase III with decreasing x in the ground state 
of the spin-orbital model G3D = (a) spin-orbital entanglement 
entropy <S° N dl-ll) and the spin-orbital correlation function 
Ci (14.11) ; (b) spin Si (14.21) and orbital Ti (I4.3j) correlation 
functions. The onset of phase V is signaled by a step-like 
increase of the entropy to <S° N = 1. Parameters: A = 0.5, 
y = —0.5 and L — 8. 


correlation function on a bond [19], 


Ci = lY, [((^•5 i+ i)(f r f m )}-(5 i 4+i)(f r f, +1 ) , 

i=1 

(4.1) 

and we compare it with the conventional intersite spin- 
and orbital correlation functions: 



The above general expressions imply averaging over the 
exact (translation invariant) ground state found from 
Lanczos diagonalization of a ring of length L. While 
Si (j4.2p and T\ (j4.3p correlations indicate the tendency 
towards particular spin and orbital order, C\ (ED quan¬ 
tifies the SOE — if Ci 7 ^ 0 spin and orbital degrees of 
freedom are entangled and the mean-field decoupling can¬ 
not be applied in Eq. ED as it generates systematic 
errors. 

To gain a better insight into the nature of a phase 
transition between phases II (AF/FO) and V and be¬ 
tween V and III (AF/AO) which occur for decreasing x 
at a fixed y < — 1/4, we study SOE vNE <S^ N , joint spin- 
orbital Ci ED, and individual spin S\ (14.21) and orbital 
T\ £OJ) correlations in Fig. [TJ Both S® N and Ci show 
a very similar behavior with a maximum within phase 
V, see Fig. [Tfa). The SOE is lower in phase III than 
in phase V, i.e., below x ~ —0.25, and decreases further 
with decreasing x. These two phases have rather similar 
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spin correlations Si, but orbital correlations T\ are sim¬ 
ilar to spin ones only within phase III (AF/AO); above 
x ~ —0.25 they vary fast within phase V and almost 
disappear (Ti ~ 0) near the transition point to phase II 
(AF/FO), see Fig. Efb). 


B. Isotropic spin-orbital SU(2)®SU(2) model 

The phase diagram of the isotropic SU(2)®SU(2) 
model (Fig. [5]) includes the same six phases as the one 
obtained at A = 0.5 (Fig. [5]), with spin-orbital corre¬ 
lations explained in Table I. The main differences to the 
phase diagram at A = 0.5 is a somewhat reduced stability 
range of phase IV (FM/AO), and also phase II (AF/FO), 
both destabilized by enhanced spin-orbital fluctuations in 
and around phase III (AF/AO). The range of entangled 
ground states is broad and includes phase III as well as 
phases V and VI on its both sides. The phase transitions 
from phase I where all quantum fluctuations are absent 
to either phase II or IV are given by straight lines and 
may be determined using mean-field approach. 

It is quite unexpected that the dimerised phases V 
and VI survive in the phase diagram of the isotropic 
SU(2)®SU(2) model at A = 1.0, see Fig. [51 These two 
phases emerge in between a disentangled phase II and 

III in the case of phase V, and similarly between phases 

IV and III in the case of phase VI, and are stabilized by 
robust quadrupling of orbital or spin correlations which 
was overlooked before [43j. This is not so surprising as 
one expects that isotropic spin-orbital interactions would 
lead instead to uniform phases only. In each case the ef¬ 
fective exchange interaction changes sign in one (either 
spin or orbital) channel which resembles the mechanism 
of exotic magnetic order found in the Kugel-Khomskii 
model [85| . 

The phases V and VI emerge by the same mechanism 
as phase V for the Ising orbital interactions, see Sec. IIII1 
In the isotropic model this phase and phase VI with com¬ 
plementary spin-orbital correlations occur in a symmetric 
way with respect to the x = y line, see Fig. [51 Orbital 
correlations in the case of phase V (spin correlations in 
the case of phase VI) change gradually towards FO (FM) 
order in phase II (IV) with increasing x ( y ). The qua¬ 
drupling of the unit cell seen in both phases in such cor¬ 
relations, shown in Fig. [9j may be seen as a precursor of 
this transition. Both phases are stable only in a rather 
narrow range and disappear for y — oo or x —» — oo, 
respectively, as presented in the inset of Fig. [5] for phase 
VI. The spin-orbital interactions and the mechanism sta¬ 
bilizing these phases are different from spin-Peierls and 
orbital-Peierls mechanisms in the ID SU(2)®SU(2) spin- 
orbital model with positive exchange (J = —1) [63|. We 
emphasize that the present mechanism of dimerization is 
effective only in one (spin or orbital) channel and thus 
it is also distinct from the spin-orbital dimerization in a 
FM chain found at finite temperature [86j |. 

A characteristic feature of SOE in phase VI (and simi- 
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FIG. 8. (Color online) Phase diagram in the (x,y) plane and 
spin-orbital entanglement <S° N (jl-lD (right scale) in the ground 
state of the isotropic spin-orbital SU(2)®SU(2) model (12.11) 
with A = 1.0 obtained for a ring of L = 12 sites. Phases 
I-IV correspond to FM/FO, AF/FO, AF/AO, FM/AO order 
in spin-orbital sectors. Orbitals (spins) follow a quadrupled 
pattern accompanied by spin (orbital) dimer correlations in 
phase V (VI). The phase boundaries determined by dominant 
modes of structure factors, shown by solid (dashed) lines, are 
of first (second) order. Inset shows the extended range of 
phase VI which separates phases IV and III for —10 < x < 0. 


lar in phase V) is a competition between the spin (orbital) 
quadrupling correlations along the chain which support 
orbital (spin) singlets, with the AF/AO order character- 



FIG. 9. (Color online) The structure factors obtained for a 
spin-orbital ring Eq. (12.11) of L — 8 sites in the full Hilbert 
space at A = 1.0: (a) spin S zz (k ), and (b) orbital T zz (k). 
The points (0.6,—1.0) (filled squares) and (—1.0, 0.6) (filled 
circles) correspond to phase V and VI, see Fig. [8j 
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FIG. 10. (Color online) Spin-orbital entanglement entropy 
<S°n in the ground state of the SU(2)0SU(2) spin-orbital 
model <EU as a function of y for representative values of 
x £ [—2.0,—0.3]. The onset of phase VI at decreasing y is 
detected by a step-like increase of entropy at the IV-VI phase 
transition. 



FIG. 11. (Color online) von Neumann entropy S° N (H-1D and 
joint spin-orbital bond correlation C\ ED as obtained in the 
SU(2)(g)SU(2) model ( 12 . 11 ) along the symmetry x — y line 
in the phase diagram with the ring of L = 8 sites for: (a) 
the model with negative coupling — J — — 1 1436 ] . and (b) the 
model with positive coupling —J — 1 constant [ 561 ] . 


istic of phase III (AF/AO). For x £ [—1.0, —0.3] the vNE 
increases discontinuously at the phase transition IV- 
VI, next drops somewhat and next increases further, see 
Fig. [TOj It exhibits a broad maximum, moving to lower 
values of y with decreasing x. This behavior shows that 
two phases (VI and III) compete in this regime. For still 
lower values of x the SOE entropy is smaller and almost 
constant when y decreases deeply into phase III. 


C. Entanglement in the SU(2)®SU(2) models 

The ground state obtained in the present spin-orbital 
SU(2)®SU(2) model for phase III (AF/AO) is distinct 
from the one found for positive coupling constant, i.e., 
J = — 1 in Eq. (12.ip . We elucidate this difference by 
studying both models along the symmetry line x = y 
in the phase diagram. Phase I (FM/FO) is found in 
the present case for x = y > 1/4, while in the case of 
positive coupling constant it becomes the ground state 
for x = y < —1/4 [6lj |. 

First we consider SOE detected by the vNE S^ N Eq. 
(ED, and by the joint spin-orbital correlation function 
Cl ED, see Fig. [TlJ To understand better the tran¬ 
sition from phase I to phase III, we consider the corre¬ 
lation functions along the symmetry line x = y. Phase 
I with FM/FO order is disentangled in both cases. At 
the quantum phase transition to phase III, signalled by a 
rapid increase of both S° N and 81 C\ |, we observe that the 
entropy reaches the highest value at the onset of phase 
III, and then decreases when x decreases and one moves 
deeper into the entangled phase III. In the present model 
m one finds C\ > 0 which is imposed by the negative 


coupling constant. The entropy maximum and also the 
maximum of 8Ci are sharp indeed and signal the onset 
of phase III, see Fig. da). The model with a pos¬ 
itive coupling constant behaves differently — here the 
joint spin-orbital correlations are negative C\ < 0, and 
both S° N and 8|Ci| have flat maxima in a range of x 
and only at x ~ 0.5 both drop rapidly. This large SOE 
for x £ [—0.25,0.5] indicates a spin-orbital liquid phase 
which forms near the SU(4) point x = y = 1/4 [6l|. Only 
at x > 0.5 the strong spin-orbital fluctuations are weak¬ 
ened when phase III is approached and SOE decreases. 

A special feature of the present SU(2)®SU(2) spin- 
orbital model is a very distinct behavior along the I- 
III phase transition line x + y = 1/2, see Fig. [8j The 
ground state energy of phase I (FM/FO) in which quan¬ 
tum fluctuations are absent, E° = — J( 1/4 +x) 2 , is found 
by taking exact classical values of spin dm and or¬ 
bital (Hd) correlations on the bonds, Si = T\ — 1/4. 
The energy decreases when the transition at x = 1/4 is 
approached. On the contrary, coming from the other 
side it is not allowed to assume classical correlations, 
Si = T\ = —1/4, as then the Hamiltonian would van¬ 
ish at the transition. In fact the energy Eq = — J/4 can 
be also obtained mainly from enhanced joint spin-orbital 
correlations, ((S j -S j + 1 )(T r f j + 1 )) = 5/16 and C x = 1/4, 
see Fig. |TTJ At the phase transition the spin and orbital 
correlations are very weak , i.e., Si = T\ ~ —1/8. This 
is a very peculiar situation as joint spin-orbital correla¬ 
tions cannot be factorized and damp to a large extent 
individual spin and orbital correlations. 

A qualitatively different nature of the SOE in both 
SU(2)®SU(2) models is also captured by its effect on the 
individual spin (orbital) correlations, see Fig. [12j In the 
present model (12.ip with J = 1, C\ > 0 damps individ- 
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FIG. 12. (Color online) Spin Si (14.2D and orbital T\ (14.3[) 
(Ti = Si) bond correlations as obtained in the SU(2)®SU(2) 
model (12.ID along the symmetry line x = y in the phase di¬ 
agram with the ring of L = 8 sites for: (a) the model with 
negative exchange J — 1 [43], and (b) the model with positive 
exchange J — — 1 [56[. 


ual spin and orbital fluctuations near the quantum phase 
transition stronger than in the case of positive coupling 
(J = —1), cf. Figs. fl2T a) and fl2T b). However, the joint 
fluctuations C\ decrease fast when J = 1, while they 
are robust in the spin liquid phase for J = —1 when 
x G [0.25,0.50). At x = 0 one finds S± = Ti ~ —0.37, 
already much below the value of Si = T\ ~ — 0.22 found 
for J = —1. The spin Si and orbital T\ correlations be¬ 
come similar in both phases deeply within phase III, as 
found by comparing these values at x = — 1 for J = 1 



L 

FIG. 13. (Color online) von Neumann entropy <S° N (11-ID ob¬ 
tained in the isotropic SU(2)®SU(2) model (12.ID at the max¬ 
imum seen in Fig. EH a) at x = y = 0.249 (squares) and at 
J = — 1 ai x = y = —0.249 (circles), i.e., at the onset of phase 
III, both for rings of L = 4, 6, • • • ,12 sites and the linear fits 
to the data (dashed lines). 


and at x = 1 for J = — 1. Here SOE is weak because 
spins and orbitals fluctuate almost independently. 

The sharp peak of <S^ N found near the phase transi¬ 
tion for J = 1 is rather unusual, see Fig. flll a). In 
this regime of parameters the ground state energy Eo 
is lowered by positive joint spin-orbital correlations Ci, 
while negative Si = Ti ~ —1/8 [Fig. fl2T a)] increase it 
somewhat. Large vNE S® N is found near the phase tran¬ 
sition for rings of even length, starting from ~ 1 for 
L = 4. It is remarkable that the vNE at the maximum 
scales with system size L, see Fig. [T3j This behavior is 
unique and proves that SOE which takes place at every 
bond is extensive and extends here over the entire ring. 
The model with positive coupling constant, J = — 1, has 
also a similar linear scaling, <S{J N oc L, but the entropy is 
smaller and systematic fluctuations between the rings of 
length of 4n and 4n + 2 sites, seen in Fig. [13l are distinct 
and indicate a crucial role played here by global SU(4) 
singlets. 

V. ENTANGLED ELEMENTARY 
EXCITATIONS: FM/FO ORDER 

A. Analytic approach 

When the ground state is disentangled, SOE is gen¬ 
erated locally in excited states (g?) and would not scale 
linearly with system size. Indeed, we analyzed the low- 
energy excitations of the disentangled FM/FO phase in 
the ID SU(2)®SU(2) spin-orbital model in Ref. and 
found a much weaker dependence on system size. We 
investigated the SOE for spin-orbital bound states (BSs) 
and spin-orbital exciton (SOEX) state and found a log¬ 
arithmic scaling, while the entropy saturates for other 
separable (trivial) spin-orbital excitations. One finds 
that the vNE is controlled by the spin-orbital correla¬ 
tion length £ and decays logarithmically with ring length 
L. 

Here we consider again the FM/FO disentangled 
ground state |0), obtained for the anisotropic spin-orbital 
SU(2 )®XXZ model with the exactly known ground state 
energy E°, 

H\Q) = E°\0). (5.1) 

Using equation of motion method one finds spin 
(magnon) excitations with dispersion 

us(Q) ~ (q + y) 4 - cos *2)> ( 5 - 2 ) 

and orbital (orbiton) excitations [HI, 

wt(Q) = Q + ^ (1 - AcosQ). (5.3) 

The spin-orbital continuum is given by 

Q.(Q, q) - uj s (j - q) + lo t + q) ■ 


(5.4) 
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Next, we consider the propagation of a magnon-orbiton 
pair excitation along the FM/FO chain, by exciting si¬ 
multaneously a single spin and a single orbital. The 
translation symmetry imposes that total momentum Q = 
2mir/L (for mm 0, • • • , L — 1) is conserved during scat¬ 
tering. The scattering of magnon and orbiton with initial 
(final) momenta — g, §■ +q} ({^ — q', % + #'}) and the 
total momentum Q is represented by the Green’s function 


G «'"> = z£ (5 - 5) 


for a combined spin {S Q and orbital (Tq + ^) excita¬ 


tion. The analytical form reads 

G(Q,w) = G°(Q,w) + II(Q,w), 

n(Q,w) = - 


2(1 + A) + (1 — A 2 )F SS (Q,to) u2 


4[1 + A (Q, w)] 


(5.6) 


2(1-A) + (l-A 2 )F cc (Q,lo) - 


+ 


4[1 + A(Q, w)] 
(1-A 2 )F sc (Q, w ) 


h;,(q,u>) 


2[1 + A(Q, w)] jffcc(Q ’ W) ’ (5J) 
where the noninteracting Green’s function is given by 

G°(Q,Lo) = jY / G° qq (Q^), (5.8) 


(Q,a;) = 


w - £l(Q,q) ’ 


(5.9) 


and 


H cc {Q,u) = 1 ^ (cos^ _cos <?) G° qq (Q,uj), (5.10) 

<7 ^ ' 

H ss (Q,u>) = 1 y ( sin § _sin ^) G° g (<3,w). (5.11) 
<7 ^ 7 

One finds the denominator in n(Q,co’) (15. 7|) . 


1 + A (Q, cj) 

! + -(! + A)F cc (Q, cj) 


- i (l-A 2 )F s 2 c (Q,c), 


1 + -(1 — A )F ss {Q,cd) 


with 


F fo .A 1 v (cos f-cos g) 2 
Fcc(Q,u) — — 2_^ w _ n(Qj9) > 


(5.12) 


(5.13) 


F rn .a - 1 v (sin f-sing) 2 

w -n ( Q, g) ’ 


(5.14) 


F sc (Q,oj) = -]T 


1 (sin ^ — sin g) (cos ^ — cosg) 


cj - fl(Q,g) 


(5.15) 


Here we define a phase S G [0, 27r] which quantifies the 
difference of dynamic properties of magnon and orbiton 
excitations throughout the Brillouin zone in the form 


tan 5 = 


(4 y + 1) — {Ax + 1)A 
{Ay + 1) + {Ax + 1)A 


tan 



(5.16) 


For the symmetry line x — y Eq. (|5.16)) greatly simplifies 
for the isotropic model at A = 1 and gives S(Q) = 0, 
which has been studied in Ref. [43[. Also, a quantity 
describes the position relative to the continuum is given 

by 


a{Q,u) 


lu — {x + y + \) 

HQ) 


(5.17) 


with 



The excitations at the top and at the bottom of the con¬ 
tinuum correspond to a{Q,u) = 1 and —1, respectively. 
In the noninteracting case, we have the imaginary and 
real parts of Eq. (EH), 

9G°(Q,w) = 
sRG 0 (Q,w) = 


where 6(x) is Heaviside step function whose value is zero 
for negative argument and 1 for nonnegative argument. 
The frequency dependence of the imaginary part exhibits 
square-root singularities in Eq. (I5.19p at the bottom and 
at the top of the continuum [90j . 


(9(1 - |a(Q,cj)|) 


KQW i- a 2 (Q,w) ’ 

9(a(Q,u>) - 1 ) 
b(Q)y/a*(Q,u)=l 
0 {-1 - a(Q,ui)) 


(5.19) 


(5.20) 


B. Numerical studies 

We note that the inclusion of the spin-orbital attrac¬ 
tion will smear out the singularities by Eq. (ETBl) since a 
more pronounced divergence of the numerator than the 
denominator occurs when a(Q, cd) ±1, and they cancel 
each other, i.e., ^sG{Q,uj) = 0. Furthermore, the poles 
of G{Q,uj) are determined by 

1 + A(Q, id) = 0. (5.21) 

Our analysis shows that for given Q most the real so¬ 
lutions of Eq. (15.21)) are interspersed within the contin¬ 
uum, but these modes are unstable to two free waves. 
However, a small number of solutions may lie well below 
the continuum. To investigate the spectra we begin with 
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the asymmetric SU(2 )®XXZ model. As it is shown in 
Fig. [TT] the attractive interactions shift spin-orbital BSs 
outside the continuum (9ll-[93j. The binding energy ap¬ 
proaches zero for the isotropic SU(2)®SU(2) model, but 
is finite for anisotropic SU(2 )<g>XXZ model due to a gap 
in the orbital excitation spectrum, and the small-g be¬ 
havior of the binding energy reveals that the BSs appear 
for arbitrarily small wave number. 

The BSs can be also obtained by the equation 
of motion method for a spin-orbital joint excitation, 
0). The collective mode follows from Eq. (|5.21[) . 
Such a collective spin-orbital excitation (bound state) in¬ 
volves spin and orbital flips at many sites and can be 
written as follows, 

l*(Q)> = T^a^Oe^-T-^lO) 

V rn,l 

= EV?i _7 q + J°)- (5.22) 

q 

with the coefficients 

a g = -E^ aK Q) e -(#- 9 )*. (5.23) 


The correlation length £ = l\ai\ 2 defines the average 
size of spin-orbital BSs or excitons and is much smaller 
than the system size, i.e., 0 < £ « L. The correlation 
length becomes extensive for a trivial continuum state, as 
shown in Fig. [15] The analytic solution of this equation 
is tedious but straightforward. The dispersion of the col¬ 
lective excitation, cjbs(Q )5 can be analyzed in a simple 
way at some special points, including Q = 0 and Q = tt. 

In the isotropic model (at A = 1), Eq. (|5.21|) reduces 
to 1 + F cc (Q,uj) = 0. In this case there is at most one 
solution for every Q [43|. Nevertheless, the anisotropic 
orbital coupling will induce more branches in part of Bril- 
louin zone (see Fig. [14]). When A < 1 and Q = 0, 
F sc (0,u;)=0 and then Eq. (|5.21|) reduces to 

1 + 1(1 + A)F cc ( 0, w) «= 0, (5.24) 

or 

l + l(l-A)F„(0,w) = 0. (5.25) 

The solution that follows from Eq. (|5.24 | ) is given by 


<^BS,l(0) 


_ (a + 3/3) 2 - y/(a- 17/3)(a - fi) 3 
8/3 

+ x(l-A)-3/3+i (5.26) 


where a = ( Ax + y), /? = (1 + A)/4. The instability of 
such a mode given by cjbs,i(0) = 0 sets up the threshold 
of the FM/FO state separating it from the AF/AO state 
(phase III) [H,[92|. Moreover, the solution for Eq. (|5.25|) 
is given by 

(ce + /3) 2 

^ bs , 2 (0 )=xFy~Y^Fp. ( 5 . 27 ) 



Q/k 



Q/tt 


FIG. 14. (Color online) (a) Excitation spectra of a ring of L = 
40 sites for the spin-orbital SU(2 )®XXZ model with A = 0.5 
as function of momentum Q at: (a) x = y = 1/4, (b) x = 
0.375, y = 0.125, and (c) x = 0.27, y = 0.23. The dotted 
(blue), dashed (green) lines inside the spin-orbital continuum 
Q(Q,q) denote the orbital and SOEX excitations, i.e., ujt(Q) 
and Cg>soex(Q), respectively, that are degenerate. The (red) 
solid lines show spin excitations. 


We find that when a < (1 — 3A)/4, both cjbs,i(0) and 
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FIG. 15. (Color online) The correlation length £ for increasing 
momentum Q as obtained for a ring of L = 80 sites in the 
anisotropic SU(2 )<g>XXZ model with A = 0.5. 


^bs ,2 (0) exist, while only cjbs,i(0) survives when (1 + 
A)/4 > a > (1 — 3A)/4. Finally, for a > (1 + A)/4 no 
bound states exist. 

When A < 1 and Q = 7r, one finds that F sc (tt : cj)= 0. 
Analogously, Eq. (|5.21|) has two solutions, u bs,i(tt) and 
cjbs, 2 (tt) when — (3 + A)/4 < y — Ax < (1 — A)/4, with 
explicit expressions for their energies: 


/ \ 1 1T A 

^bs,i(tt) = X + y + --— 

_ [(y+l) - A ( g +I )] 2 

1 +A 

, , 1 

wbs, 2 (tt) = X 4 - y + - 


c - (27 - 1 - A) § V27 - 9 + 9 A 
+ 8(A = 1) 


(5.28) 


(5.29) 


with 


FIG. 16. (Color online) Excitation spectrum for the spin- 
orbital SU(2)0Z2 model (|2. ID on a ring of L = 40 sites. The 
dotted (blue) and dashed-dotted (gray) lines indicate the or¬ 
bital excitation cjt(Q) and the lower boundary of the con¬ 
tinuum u c (Q) — both are degenerate. The dashed 

(green) and solid (red) lines denote the SOEX excitations 
^soex(Q), and spin excitations ojs(Q)- Parameters: A = 0.0, 
x — 0.375, and y = 0.125. 


[HI, e.g. S^Tf |0) at site /, within the spin-orbital con¬ 
tinuum with £ = 0 (see Eq. (|5.23|0 . It is straightforward 
to derive, 


[H,SfT{-] |0> = [Ci(x,y) + D t (x,y)] |0), (5.31) 


where 


Ci(x,y) — (x + y)S t T { (S l _ 1 T l _ 1 + S l+l T l+1 ) 


Di(x,y) = 


A ( x ~\) ( s i T i ' ' 5 ' T m) 




(5.32) 


7 = (y + 1/4) - A(x + 1/4), 

( = 3 + 47-47 2 - 6A-47A + 3A 2 . (5.30) 

When y — Ax < — (3 + A)/4 and (1 — A)/4 < y — Ax < 
(3A + l)/4, only one cjbs,i(tt) solution exists. In case of 
y—Ax = (1 —A)/4, cjbs ,2 (tt) merges with lower boundary 
of the continuum. 


C. Propagating spin-orbital exciton states 

Especially at the SU(4) symmetric point, i.e., at x = 
y = 1/4, spinon and orbiton are strongly coupled to 
form a joint SOEX state inside the spin-orbital contin¬ 
uum across the whole Brillouin zone. Usually such an 
elementary spin-orbital excitation in the continuum is 
unstable and decays into a spinon and an orbiton [38| . 
However, it is surprising that such a SOEX state propa¬ 
gates here as a undamped on-site spin-orbital excitation 


The dissipative term Di(x , y) vanishes when x = y = 1/4. 
Consequently, the dispersion of the SOEX state is given 
by 

<^soex(Q) = - (1 — A cos Q ). (5.33) 

The SOEX state for x = 1/4 is degenerate with the 
orbital wave excitation, see Eq. (|5.3p and Fig. m a). 
When A = 1 and x — y > 0.25, there is a quasi-SOEX 
state inside the spin-orbital continuum for Q < 7r, with 
£ < 1. When x j- y, the residual signal of the SOEX state 
denoted by finite £ vanishes and the BS appears. 

The smaller the A is, the more values of Q give cjbs, 2 (Q)- 
There is no BS at Q = 7r for y — Ax > (3A + 1)/4. 
Furthermore, away from the symmetric point, the SOEX 
will acquire a finite linewidth due to residual interactions 
into magnon-orbiton pairs, 

r = Q{G- 1 (Q, w )}. 


(5.34) 
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The exciton spectral weight can be calculated from the 
self-energy E(Q,u;), 


Z Q = 


1 - 


/ dE(Q,u) 
V foo 


-l 


(5.35) 


^soex(Q) 

where the SOEX energy is given by by the pole, 

^soex(Q) = 5R{G _ 1 (Q,u;)}. (5.36) 


If T /[oj — cjsoex(Q)] 2 —>• 0, the exciton is stable. The 
decay rate of the SOEX increases with growing x > 1/4 
and also for decreasing momenta Q, and they coincide at 
A = 0. 

In the limit of Ising orbital interactions (A = 0), the 
orbital part becomes classical and orbitons are disper¬ 
sionless, indicating localized orbital excitations, see Fig. 
m We find that two BSs, cjbs,i( 0) and cjbs, 2 ( 0 ), exist 
when y < 1/4, and there are no BSs otherwise. In con¬ 
trast, at Q = 7 r one finds two solutions, cjbspM and 
cjbs ,2 (tt) 5 when —1/4 < y < 1/4 and no BS is found 
for y > 1/4. In Fig. [16] both BSs, cjbs,i and cjbs, 2 , 
are undamped. In this case, ujt(Q) = 1/4 + x and is 
degenerate with the lower boundary of the continuum 
cj c (Q) = 1/4 + x. Moreover, cjsoex(Q) = x + y, and 
especially when x = y = 1/4, these excitations coincide. 
In this case all excitations are dispersionless, see Fig. [16] 
The spin-orbital BS in Fig. [16] appears below the bot¬ 
tom of the spin-orbital continuum and is stabilized by its 
binding energy. 


VI. VON NEUMANN ENTROPY SPECTRA 

A. The spectra in the FM/FO phase 

To investigate the degree of entanglement of spin- 
orbital excited states, we introduce the vNE spectral 
function in the Lehmann representation, 

<S v n(Q,w) = - X Tr {/ 9 s l) l°g 2 - u n (Q)} , 

n 

( 6 . 1 ) 

where we use a short-hand notation (/i) = (Q,uj n ) for 
momentum Q and excitation energy c j, and 

p { s n) =Tr T \* n (Q))(* n (Q)\ (6.2) 




FIG. 17. (Color online) The vNE spectral function <S v n(Q,Cl;) 
<nnj for the anisotropic spin-orbital SU(2 )®XXZ model 
(12.11) with A = 0.5 on a ring of L — 160 sites, see Fig. fl4l bb 
obtained for different momenta Q < 0.87T, and for: (a) even 
excitations, and (b) odd excitations. Isolated vertical lines 
below the continuum indicate the BS, with dispersion given 
by the dashed (red) line. Parameters: x — 0.375, y — 0.125. 


We have derived an asymptotic form of the vNE as a 
function of £ [43^ . 



One finds the asymptotic logarithmic scaling of vNE of 
spin-orbital BSs and SOEX state whose magnon-orbitons 
correlation length is short-range, and the vNE is given by 


is the spin (S) density matrix obtained by tracing over 
the orbital (T) degrees of freedom. In Fig. [l7]we present 
the analytic results for the vNE spectral function when 
A = 0.5, x = 0.375, and y = 0.125. The parity symme¬ 
try is broken at x 7 - y and A / 1 . The even and odd 
excitations show diverse behavior of their entanglement, 
as is displayed in Fig. [TT] Inspection of the vNE spectra 
shows that the entanglement reaches a local maximum 
at the BSs and SOEX states, and all these states have 
short-range correlation length £. 


<S vN = log 2 L + c 0 . (6.4) 

In particular, Co = 0 for the SOEX state and Co < 0 
otherwise. Such relation is displayed in Fig. [18] Fig. 
H5l implies that cjbs,i are always stable for all momenta 
while cjbs ,2 are undamped for large momenta. 

When both x and y are away from 1/4, the SOEX state 
is unstable and decays into a spinon and an orbiton. In 
this case the correlation length £ becomes extensive. We 
have verified that the scaling is entirely different in such 
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FIG. 18. (Color online) Scaling behavior of the entanglement 
entropy <S v n of the spin-orbital BSs at Q — 0.87T (points) 
for x — 0.375, y — 0.125, and for: (a) A = 0.5 and (b) 
A = 0. Lines represent logarithmic fits to Eq. (16.41) with: (a) 
co = -0.421 and -0.842; (b) c 0 = -0.368 and -1.122. 

a case and the entropy of the SOEX scales instead as a 
power law, 

<S V N = + C °- (6-5) 

The vNE saturates in the thermodynamic limit. 


B. RIXS spectral functions in the FM/FO state 

The entanglement spectral function <S v n(Q,^) has a 
similar form as any other dynamical spin or charge cor¬ 
relation function. There is, however, an important dif¬ 
ference — as there is no direct probe for the vNE of an 
arbitrary state, the SOE spectra can be calculated but 
cannot be measured directly. On the other hand, we 
have shown before [43] that the intensity distribution of 
certain RIXS spectra of spin-orbital excitations in fact 
probe qualitatively SOE. 


We introduce the spectral function of the coupled spin- 
orbital excitations at distance /, 

MQ,cj) = - lim Im(0|rg )t , - 1 - . rg } |o). 

7 T ?7-^0 ^ UJ + Eo — H — IT] ^ 

( 6 . 6 ) 

Here 

T Q=jlY,e ' 0 ’ S ; T , 7 < 6 - 7 > 

is the local operator for an on-site joint spin-orbital exci¬ 
tation measured in RIXS [ini looj . We employ as well the 
even and odd parity operators, 

r r ) =^E e * i ( s « ±s H) i ;‘. (“> 

which serve to probe the nearest neighbor spin-orbital 
excitations. 

Intuitively, the on-site spectral function Ao(Q, uj) high¬ 
lights the SOEX state. It is found that both cjsoe x(Q) 
and cubs(Q) are solutions of Eq. (|5.21|) when x = y = 
1/4. However, the weight of the BS is here zero. Now 
the spectral function is given by 

A 0 (Q,lo) = S {uj - cjsoex(Q)} , (6.9) 

which is confirmed in Fig. [T9l a). As the point (x,y) 
moves away from the symmetric point, i.e., x = y = 1/4, 
the BSs gain spectral weight which decreases with mo¬ 
mentum Q. The spectral weight vanishes at Q = tt ac¬ 
companying a square-root singularity at the lower bound 
of the continuum, see Fig. [T9l b). The S peak turns into 



FIG. 19. (Color online) The spectral function of the 
on-site spin-orbital excitation Ao(Q,cj ) for the anisotropic 
SU(2 )<g>XXZ model with A = 0.5, and for: (a) x = 0.25, 
y — 0.25, (b) x — 0.5, y — 0.5, and (c) x — 0.375, y — 0.125. 
In each panel the momenta Q range from 7r/10 (bottom) to 
97 t/ 10 (top); the peak broadening is rj = 0.001. The red dot¬ 
ted (green dashed) fines mark the positions of BSs (SOEX) 
states. Gray dash-dot fine signals the onset of the continuum. 
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FIG. 20. (Color online) The spectral function of the spin- 
orbital excitation at nearest neighbors Ai+(Q,lj) with even 
parity for the anisotropic SU(2 )<g>XXZ model with A = 0.5, 
and for: (a) x — 0.25, y — 0.25, (b) x — 0.5, y = 0.5, and (c) 
x = 0.375, y = 0.125. In each panel the momenta Q range 
from 7 r /10 (bottom) to 97 r /10 (top); the peak broadening is 
77 = 0.001. The red dotted (green dashed) lines mark the 
positions of BSs (SOEX) states. Gray dash-dot line signals 
the onset of the continuum. 


a broad peak. A difference between x and y induces a 
second branch of BSs and it gains larger spectral weight 
at large momenta than the first BS, see Fig. HW c). Alto¬ 
gether, the evolution of spectral weight is similar to that 
of the correlation length in Fig. [15l 

The BSs can be captured also by the spectral func- 



00 CO CO 


FIG. 21 . (Color online) The spectral function of the spin- 
orbital excitation at nearest neighbors Ai+(Q,u) with odd 
parity for the anisotropic SU(2 )®XXZ model with A = 0.5, 
and for: (a) x — 0.25, y — 0.25, (b) x — 0.5, y — 0.5, and (c) 
x — 0.375, y — 0.125. In each panel the momenta Q range 
from 7r/10 (bottom) to 97 t/ 10 (top); the peak broadening is 
77 = 0.001. The red dotted (green dashed) lines mark the 
positions of BSs (SOEX) states. Gray dash-dot line signals 
the onset of the continuum. 


tions for the nearest neighbor excitations, see Figs. l20l 
and [2TJ With increasing momentum Q, the spectral 
weight of even-parity excitation Ai+(Q,u;) at the first 
BS decreases (see Fig. [20|), while the spectral weight 
of odd-parity excitation Ai-(Q,lj) rises (see Fig. |2T|) . 
Ai_|_(Q,a;) reaches its valley for quasi-SOEX states, but 
no such feature could be found in Ai_(Q,cj). 


VII. DISCUSSION AND CONCLUSIONS 

Motivated by the discovery of new Majumdar-Ghosh- 
like valence-bond spin singlet phases triggered by or¬ 
bital correlations, we have studied the spin-orbital entan¬ 
glement (SOE) in the one-dimensional (ID) anisotropic 
SU(2 )®XXZ spin-orbital model with the negative ex¬ 
change interaction. The asymmetry between spin and 
orbital degrees of freedom yields a better insight into the 
phase diagram and the mechanisms responsible for the 
different types of order observed for this system. In ad¬ 
dition to the four uniform phases I-IV, our study demon¬ 
strates that a gapful phase V exists in case of classical 
Ising orbital interactions, i.e., in the SU( 2 )(g)Z 2 model. 
It is characterized by quadrupling of the unit cell seen as 
a maximum of the orbital structure factor at k = tt/2. 
For y = —1/4 this provides a perfect dimer structure 
of spin singlets in the whole region of stability of this 
phase, where the dimer spin correlations D(r) develop 
and uncover long-range dimer order. The dimer phase V 
is quite robust and survives when the orbital quantum 
fluctuations at A > 0 are taken into account. 

The phase diagram is still richer at finite A > 0, when 
quantum orbital fluctuations develop and induce an or¬ 
bital dimer phase VI, with a complementary role of spin 
and orbital correlations to phase V. The emergence of 
the nonuniform phase V is a result of the joint interac¬ 
tion between spin fluctuation and orbital degree of free¬ 
dom, and thus phase V carries finite SOE. The orbital 
fluctuations enhance the SOE in ground state near the 
III/V phase transition and lead to phase VI when the 
{x,y} parameters are interchanged. We also anticipate 
that these dimer phases may survive in higher dimen¬ 
sions. In fact, Lieb-Schultz-Mattis theorem is applicable 
to the present model for arbitrary x and y and nonzero A, 
where a finite gap exists above these degenerate ground 
states. Both phases V and VI are gapped phases with 
alternating spin and orbital singlets, respectively. As we 
have shown, the phase boundaries can be captured by 
SOE and fidelity susceptibility. The phase transition be¬ 
tween phases III and V is a first-order transition in the 
A = 0 case, and the transition changes to continuous 
when A > 0. 

An important consequence of finite SOE in the ground 
state is that it invalidates mean-field decoupling of spin 
and orbital degrees of freedom, as this would imply a 
spin-orbital product ground state. A similar restriction 
applies to the entangled elementary excitations in the 
disentangled ferromagnetic phase with ferro-orbital or- 
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der in spin-orbital systems which were analyzed he: 
with the help of the von Neumann entropy spectral fun 
tion. Spin-orbital excitations are highlighted by nontrr 
ial SOE, especially by logarithmic scaling of SOE in th 
phase. 

A priori , the SOE makes it necessary to treat the eige] 
states of a given model exactly. In fact, since a mean-fie' 
decoupling shields the SOE, it fails to describe eleme] 
tary excitations even qualitatively correctly in a min 
ber of spin-orbital models. In antiferromagnetic grour 
states with ferro-orbital order it was demonstrated bot 
in theory (3?1 and experiment [H| that the spin-orbit 
excitation fractionalizes into freely propagating spine 
and orbiton, giving rise to spin-orbital separation und< 
specific condition. The SOE in the spin-orbital separ, 
tion remains unclear. The low-lying excitations in phas< 
II (AF/FO) and III (AF/AO) are spin waves with vanis] 
ing SOE, corresponding to a two-spinon continuum of c 
antiferromagnetic spin chain. The low-lying excitation : 
phases V and VI corresponds to spin-orbital excitatio: 
as shown in Figs. [2jt>) and [2fc). The problem of tl 
SOE in elementary excitations in other phases remains 
challenge for future studies. 

Summarizing, we have shown that the anisotrop 
SU(2 )(g)XXZ spin-orbital model with negative exchan^ 
coupling has remarkably different behavior and phase d 
agrams from the well known SU(2)0SU(2) model wil 
positive exchange coupling. While the spin-orbital li< 
uid phase is absent in the former case, we have four 
that the joint ferromagnetic/ferro-orbital fluctuatioi 
are surprisingly strong at the quantum phase tram 
tion to the antiferromagnetic spin order which gives ew 
stronger SOE than that established for the ID isotrop 
SU(2)(g)SU(2)model with positive exchange coupling. 
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Appendix: Phase diagrams for the two-site model 


To understand better the phase boundaries in the 
anisotropic SU(2 )<g>XXZ spin-orbital model with neg¬ 
ative exchange coupling dm we present here an exact 
solution for the system of L = 2 sites 0, 


%12 =~J 


(& • S 2 + x) 


- (T+T 2 - + T-T+) + TfTf + y 


(A.1) 




X 

FIG. 22. (Color online) (a) The phase diagram of 2-site 
model when A = 0. (b) A = 0.5. Phases I-IV correspond 
to FS/FO,AS/FO, AS/AO, FS/AO configurations in spin- 
orbital sectors. 


Again, {x, y} are the parameters, and 0 < A < 1 inter¬ 
polates between the Ising Z 2 (A = 0) and Heisenberg 
SU(2) (A = 1) symmetry. The orbital interaction with 
XXZ symmetry can be exactly diagonalized and one 

finds 4 eigenstates: | + +), |-), (| H—> + | — +)) /y/2, 

(| H—) — | —b))/V2, corresponding to eigenvalues: 1/4, 
1/4, A/2 — 1/4, —A/2 — 1/4, respectively. At A = 0 we 
recover doubly degenerate configurations from the latter 
two, while at A = 1 we recover a triplet T = 1 from 
the first three. In any case, the third component of the 
triplet, (| H—) + | —b)) /V% is always an entangled ex¬ 
cited state with the present choice of parameters, while 
the orbital singlet, (| H—) — | — b)) /\/2, is an entangled 
ground state for some parameters. 

The spin part is classified as a triplet S = 1 or a singlet 
S = 0, and these states are accompanied by the orbital 
states described above. This gives the states I-IV in Table 
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TABLE II. The spin-orbital configuration for the model (1A.1D . 
with phases I-IV defined by distinct spin S and orbital state 
as as obtained at 0 < A < 1. The lowest energy E° / J has 
degeneracy d and becomes the ground state at the respective 
values of {x, y} parameters, Fig. [22] At A = 0 the degenera¬ 
cies of phases III and IV change to 2 and 6 for the Z 2 0 Z 2 
symmetry, while at A = 1 the degeneracies for the ground 
states I-IV are 9, 3, 1, 3 and follow from the SU(2)0SU(2) 
symmetry. 


phase 

S 

orbital states 

E°/J 

d 

I 

1 

+ +)>-) 

- (3 + x ) ii + y) 

6 

II 

0 

+ +)? 1-) 

-(-§ + *) (i + y) 

2 

III 

0 

+ 

1 

1 

Hy 

H + d(D#-0 

1 

IV 

1 72 (1 + - > -1 +)) 


3 


II, and one of them is the ground state for any point in the 
(x, y ) plane. All phase transitions are first order, with a 
change of spin or orbital state. The phases II and IV are 
symmetric for A = 1, and the transition between I and III 
occurs along the x + y = 1/2 line [IJ. At A = 0 phase IV 
exists for y > 1/4 while phase II only for x > 3/4, see Fig. 


l22l a). This reflects the essential difference between the 
orbital configurations in the Ising limit and the quantum 
spin states. Note that spin singlet S' = 0 is an entangled 
state, but the larger Hilbert space gives no spin-orbital 
entanglement in any phase. Thus, the total energies and 
the phase diagram are easily deduced, see Fig. [22] 

From the comparison of the energies E° of phases III 
and IV (Table II) one can see that the boundary between 
III and IV, given by the straight line y = 1/4 + A/2, 
moves upwards with increasing A, see Fig. l22fb). Ac¬ 
cordingly, the phase boundary between phases I and III 
is also modified. The interplay between spins and or¬ 
bitals develops and leads to interesting consequences of 
entanglement for rings with L > 4. Remarkably, the 
trivial phase diagram found for Ising-Ising interactions 
for L = 2 sites (with II-III and IV-III transition lines at 
x = 1/4 and y = 1/4, similar to the diagrams in Fig. 
122 ]) is the same for the Z 2 0 Z 2 spin-orbital model in the 
thermodynamic limit, 

u M2 = ~ j E ( s j s Ui + x ) ( T J T i+i + y) • ( A - 2 ) 
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